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We investigate the long-time behavior of the survival probability of a tagged particle in a single-file diffusion 
in a finite interval. The boundary conditions are of two types: 1) one boundary is absorbing the second 
is reflecting, 2) both boundaries are absorbing. For each type of the boundary conditions we consider two 
types of initial conditions: a) initial number of particles N is given, b) initial concentration of particles is 
given {N is random). In all four cases the tagged-particle survival probability exhibits different asymptotic 
behavior. When the both boundaries are absorbing we also consider a case of a random interval length 
(single- file diffusion on a line with randomly distributed traps) . In the latter setting, the initial concentration 
of particles has the same effect on the asymptotic decay of the survival probability as the concentration of 
traps. 
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I. INTRODUCTION 

The diffusion of particles in narrow channels where the 
particles cannot pass each other is known as the single-file 
diffusion (SFD). Although the overall density of particles 
in such a system evolves as if there is no interactionji the 
diffusion of one specific particle (a tagged particle, or a 
tracer) is slowed down by the exclusion interactions with 
neighboring particles. Experimentally SFD is observed 
for instance in narrow biological pores, ^ inside zeoliteSj'^'^ 
during the sliding of proteins along DNA,^'^ and in many 
other either natural or artificial systems^^-il.. 

Theoretical description of the dynamics of tagged- 
particles in SFD was pioneered by HarriSfi^ He showed 
that the mean-square displacement of the tagged par- 
ticle grows with time as t^/^ (in contrast to the linear 
increase for the normal diffusion). After that many au- 
thors studied the tagged-particle dynamics under differ- 
ent conditions (for a comprehensive review see Introduc- 
tion in Ref . KLQl) . Despite the long history of the problem, 
studies devoted to SFD in finite systems appears only 
recently.—"— Most of these studies have focused on closed 
systems where the number of particles is conserved. An 
exception is Ref. [2j where the tagged-particle propagator 
for a finite interval with two absorbing boundaries was 
discussed as a particular application of the general for- 
malism. In the present paper we go beyond the previous 
studies in the followings: 1) we consider SFD in a fi- 
nite interval with absorbing boundaries, 2) instead of the 
moments of tagged-particle probability density function 
we study the long-time behavior of the tagged-particle 
survival probability. 

The absorbing boundary conditions can be used to 
model traps, and the trapping model (capturing of diffus- 
ing particles by immobile traps) represents the simplest 
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model of diffusion-limited reactions^ir— which occur in 
many real systems (examples of such systems are cited in 
the first paragraph). Up to our best knowledge, there are 
only a few studies of the effect of hard-core interaction on 
the first passage properties. For a recent study of a first 
passage time density of a tagged particle in an infinite 
system see Ref. [23. There are also earlier works focusing 
on the overall absorption rate in systems with dissimilar 
interacting particles. ^*'^^ The semi-infinite interval with 
the absorbing boundary at the origin was considered in 
Ref. [so, where the asymptotics of the tagged-particle sur- 
vival probability was discussed. 

In the present paper we consider a diffusion of identi- 
cal hard-core interacting Brownian particles in the one- 
dimensional interval of the length L. We discuss two 
types of initial conditions and two types of boundary con- 
ditions. As for the boundary conditions, we consider: 1) 
a perfectly absorbing boundary at x = and a reflecting 
boundary &t x = L, and 2) both boundaries are perfectly 
absorbing. As for the initial conditions, we discuss two 
settings: a) initially there are exactly N particles homo- 
geneously distributed within the interval; b) the initial 
concentration of particles c is given and it is homoge- 
neous within the interval. Hence, in the latter setting, 
the initial number of particles is a Poisson random vari- 
able with a mean value cL. In all four cases it is assumed 
that as soon as the particle hits an absorbing boundary, 
it is immediately removed from the system. 

The paper is organized as follows. In Sec. [TTl we con- 
sider general boundary conditions. First, for a given ini- 
tial number of particles N (Sec. Ill Ap . we present the 
probability density function (pdf) of a tagged particle's 
position (Eq. ([2])) and a survival probability of the tagged 
particle (Eq. (J4])). The pdf ([2]) represents the starting 
point of all further analyses. Second, we derive the same 
two quantities for a given initial concentration c (Sec. 
IIIBi Eqs. dZl), ([12])). Sec. llHl presents the asymptotics of 
the survival probability in the case when the left bound- 
ary is absorbing and the right boundary is reflecting (Eqs. 
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dlOl), dH]) for a given N, Eq. ^ for a given c). Sec.HVl 
contains the asymptotics of the survival probability in 
the case when both boundaries are absorbing (Eq. (|29|) 
for a given N, Eq. ([5T|) for a given c). Further, in Sec. 
IIVDI the survival probability is discussed when the in- 
terval length L is random (cf. Eq. ((33|)). The paper is 
concluded by the interpretation of the results fSecFV]). 



II. DYNAMICS OF A TAGGED PARTICLE: GENERAL 
BOUNDARY CONDITIONS 

A. Fixed initial number of particles 

The basic quantity which is assumed to be known is the 
pdf of a single noninteracting particle, f{x, t). It satisfies 
the diffusion equation 

-f{x,t)=D — fix,t) , xeiO,L), (1) 

subject to appropriate boundary conditions at a; = and 
X = L, and to the initial condition f{x,0) = l/L, x e 
(0,L). 

Assume that initially there are N particles distributed 
randomly in (0, L). Let us label the particles by numbers 
1, . . . , A^, from left to right. Due to the hard-core inter- 
action the particle order is conserved for all times. For 
a given x, x £ (0, L) the exact pdf of the tagged particle 
with number n, n = 1, . . . , N, reads 



Pn:Nix,t) 



m 



[n- iy.{N 



-Ji^.t) 



(2) 



x(F(x,t)r-i(l-^(x,t))^-": 



where F{x, t) is a probability of finding the single nonin- 
teracting particle to the left of x (either in a left trap or 
diffusing within an interval (0, x)). We specify F{x, t) in 
subsequent sections (cf. Eq. ()17|) when the right bound- 
ary is refiecting, and Eq. (P7|) when both boundaries are 
absorbing) , the results of the present section are valid for 
a general F{x,t). Pdf (HJ constitutes the starting point 
for all subsequent analyses. 

To justify Eq. ([2]) notice the following equivalence be- 
tween the dynamics of hard-core interacting particles and 
of noninteracting ones. If two neighboring hard-core par- 
ticles collide, they bounce off each other and their order 
is preserved. Equivalently, we can assume that, instead 
of the collision, the two particles pass freely (without in- 
teraction) through one another and after that we just 
exchange their labels. This is the main idea behind Eq. 
([2]). According to this idea the dynamics of interacting 
particle is generated merely by exchanging the labels of 
noninteracting particles whenever the two of them pass 
each other (for an illustration, see Fig.[T]). Therefore the 
pdf for the n-th tagged particle equals to the pdf for the n- 
th left-most particle among N noninteracting ones. The 
latter pdf is on the right-hand side of Eq. ©. To see 
this, assume that each of the noninteracting particles is 




FIG. 1. Illustration of space-time trajectories of A = 4 parti- 
cles. During the time interval (0, t), the tagged particle (bold 
line, n = 2) undergoes three elastic collisions with its neigh- 
bors. Each collision can be equivalently represented as fol- 
lows: we let two approaching particles pass freely through 
one another, and after they cross we exchange their labels. 
This maps the trajectories of interacting particles onto the 
trajectories of noninteracting ones. Namely, the trajectory of 
the tagged particle is equivalently represented by the trajec- 
tory of the second left-most particle among A = 4 noninter- 
acting particles (regardless the actual identity of the second 
left-most noninteracting particle). The mapping implies that, 
at the time t, the pdf for the tagged particle equals to the pdf 
for the second left-most particle among A = 4 noninteracting 
particles. For a general n, and A, the pdf is given in Eq. ((2)1. 



described by pdf f{x,t), then the right-hand side of Eq. 
^ (multiplied by dx) yields the probability that there 
is a particle in (x, x + dx) and, simultaneously, there are 
precisely n — 1 (N — n) particles to the left (to the right) 
of it. The factorials account for all possible sequences of 
particles to the left and to the right of the n-th particle. 
For a further discussion see Refs. [l|, [3l|, or [S^- The pdf 
^ can also be derived formally by solving a hierarchy 
of coupled diffusion equations as it was done in Ref. |30| 
for a semi-infinite interval with the absorbing boundary 
at the origin (cf. Eq. (35) in Ref. [13). 

The survival probability, that is the probability that, 
the tagged particle does not hit an absorbing boundary 
before t is given by the spatial integral 



Sn:N{t) = / dxpn:N{x,t) 



(3) 



It reads (see Appendix |X] for derivation) 



n J ^ A — n + K + 1 



X (1-F(0,i)) 



W-n+fc+l 



-(l-F(L,i)) 



n — 1 

Af-n+fc+l 



(4) 
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B. Fixed initial concentration of particles 

The finite interval with a fixed initial concentration of 
particles c can be modeled as follows. We imagine that, 
for t < 0, the interval (0,L) is a part of the whole real 
line on which the particles are distributed randomly with 
the uniform concentration c. Subsequently, at i = 0, the 
interval is decoupled from the line by the boundaries at 
a; = and x ^ L. During the diffusion (i > 0), if a 
particle hits an absorbing boundary, it leaves the interval 
and will never return. 

Initially there are on average cL particles within the 
interval (0,L). The probability that the initial number 
of particles equals N is given by the Poisson distribution 



P{N) 



icL) 



N 



m 



(5) 



For a given N, the tagged particle pdf is given by Eq. ^. 
In the present case, the pdf of the n-th tagged particle 
regardless the initial number of particles is obtained from 
pdf (HI by the summation 



Pn{x,t) = ^ Pn:N{x,t)P{N) , 



(6) 



N=n 



which yields 



Pn{x,t) ^cLf{x,t) 



[cLF{x,t)]' 
(n-1)! 



■exp[-cLF(a;,t)] (7) 



How to interpret pdf p„(a;,i)? Is it normalized to one 
at i = 0? To give the answers let us return for a while to 
the fixed A'' case. It is easy to show that 



Pn:N{x,t) = Nf{x,t) , / dxpn:Nix,0) = 1 . (8) 



N 

E 

n=l 



The first equation manifests the well known fact that the 
overall density of particles evolves as if there is no inter- 
action. The second equation is a normalization condition 
- we assume that, at t = 0, there are N particles in the 
interval hence the n-th particle is obviously also there. 
Returning to the fixed c (random N) case, the first equa- 
tion has its counterpart in 



oo 

E 



Pn{x,t) = cLf{x,t) 



(9) 



where cL is the mean number of particles initially dis- 
tributed in (0,i). The normalization condition under- 
goes more significant changes. When N is random, it is 
no longer sure that the n-th particle is initially located 
in (0,L). Therefore we have 



dxpnix, 0) = 1 - p{N < n) 



where p{N < n) stands for the probability that the initial 
number of particles is less than n. 



The probability that, at the time t, the n-th tagged 
particle is in (0, L), 



Sn(t) = / dxpnix,t) , (11) 

Jo 

reads (cf. Appendix E| 

s,,o. ' fH£<£f!!x (12, 

n — 1 -^-^ n + k k\ 

(^(i,i))"^''-(J^(0,t))"+'' 

Notice that the survival probability p^ can also be ob- 
tained from Eq. ^ as follows: 



Sn{t) = J2 Sn:N{t)PiN) 



(13) 



N^ 



The results derived in the present section hold true for 
any of the two types of boundary conditions. Let us now 
treat individual settings in a more detail. 



III. THE LEFT BOUNDARY IS REFLECTING THE 
RIGHT BOUNDARY IS ABSORBING 

A. Single noninteracting particle 

The pdf of the single noninteracting particle satis- 
fies the diffusion equation ([1]) subject to the absorbing 
boundary condition at the origin, /T(0,i) = 0, the re- 
flecting boundary condition at a; = L, d f^ / dx\x=L = 0, 
and the initial condition /T(a::, 0) = 1/i, x G (0,L). The 
solution of this initial-boundary value problem can be 
expressed as the eigenfunction expansion 



/t(->*)-37E 



(-1)* 



7ri ^^ 2fc + 1 
k=a 



cos 



(l-f)f(2.+ l) 



X exp 



[^(2^ + 1) 



Dt 



(14) 



Each term of the above series decays exponentially in 
time with a decay rate being higher for higher values of 
k. The survival probability of the single noninteracting 
particle. 



ST{t)= / AxMx,t) , 
Jo 



(10) has an asymptotic representation 



ST{t) 



exp 



-(^)- 



(15) 



(16) 
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B. Fixed A^ 



C. Fixed initial concentration of particles 



Let us choose an arbitrary fixed x, x £ (0,i). At 
a given time t, the probability that the noninter acting 
particle is found to the left of x, F{x,t), is given by the 
sum of probabilities of two disjoint events. First, the 
probability that the particle has already been trapped 
before t. Second, the probability that, at the time t, the 
particle is found within the interval (0, x). Since the right 
boundary is reflecting, we get 



F{x,t) = (1 - ST{t)) + / dx'h{x\t) . (17) 
Jo 



It follows that 



1 ~ F{0,t) ^ STit) , l-F{L,t) = 



(18) 



We introduce the two results (TTSt into Eq. ((4]) thus ob- 
taining the exact expression for the tagged particle sur- 
vival probability. We have 



(-1)'= 



(19) 



n-l 

E 

fe=0 



iV-n-|-fc + 1 



n- 1 



iSAt)) 



N-n+k+l 



Each summand of the above sum decays exponentially 
with time (cf. Eq. ([T61) '). The leading term in a long-time 
limit is that with fc = 0. It reads 



^n:N{t) 



iV — n + 1 V n 



^^ {ST{t)f''^^' ■ (20) 



Substituting ([TH)) into the above result we obtain the 
explicit form of the tagged particle asymptotic survival 
probability: 



where the prefactor reads 



/ TT \' 



-iN-n + l)(-)Dt 



^n:N 



N' 
N -n+l\n 



N-n+l 



(21) 



(22) 



Notice that the asymptotic tagged particle survival 
probability decays exponentially with time, the decay 
rate being determined by the number of particles located 
between the tagged particle and the hard wall, (N—n). In 
particular, the survival probability of the right-most par- 
ticle (n = N) decays exactly with the same rate as that of 
the single-diffusing particle (cf. Eq. ([TCI) ). In this regard 
the result ((20|) is similar to the one obtained in Ref. 30 for 
a semi-infinite interval with the absorbing boundary at 
the origin (where the single-particle survival probability 
decays as a power law t~^/^). 



For a given N, Eqs. (I^Hl) . (PT|) demonstrate an effect of 
entropic repulsion among the particles: the more parti- 
cles are located between the tagged particle and the re- 
flecting boundary the higher the decay rate of the tagged 
particle survival probability. Hence, in the present set- 
ting (when c is given but N is random), one could expect 
that the decay rate of the tagged particle survival prob- 
ability is controlled by an average number of particles 
located between the tagged particle and the reflecting 
boundary. However, this is not the case. 

According to Eq. ([T^ , the survival probability of the 
n-th particle is given by 



Slit) 



E 



Sl.NWPiN) 



(23) 



In the long-time limit, only the first term (N = n) of 
contributes significantly to the sum (cf. Eqs. (PO)) . (pij) ) 
Thus the asymptotic survival probability reads 



Slit) - cL 



icLrz!_ 

in-iy. 



'Srit) 



(24) 



Hence the long-time behavior of the survival probabil- 
ity S'^it) is dominated by those cases when the tagged 
particle is initially the right-most one iN = n). 



IV. BOTH BOUNDARIES ARE ABSORBING 

A. Single noninteracting particle 

Let us now turn to the second generic case when the 
both boundaries are absorbing. Now, the diffusion equa- 
tion ID) is supplemented by the boundary conditions 
/A(0,t) = /A(L,t) = 0, and by the initial condition 
/a(x,0) = 1/L, X E iO,L). The solution reads 



f^(x, t) = —Y ^^"(^^(^^ + ^)/-^) ^-l.{2k+l)/L]'Dt 

TT T , ^- ^ 



fe=0 



2fc 



(25) 

The survival probability of the single noninteracting 
particle, S'A(i), is given by a spatial integral of pdf 
fAix,t) over the interval (cf. Eq. (fT5|)). The asymptotics 
of SAit) reads 



SAit) 



exp 



(i 



Dt 



(26) 



B. Fixed A^ 

Since /a (2:, 0) is uniform, the probability that the non- 
interacting particle is trapped by the left boundary before 
t equals (1 — SAit)) /2. Hence the probability Fix,t) to 
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find the noninteracting particle to the left oix, x G (0, L), 
is given by (compare this result with Eq. (|17l) ) 

F{x, t)^^{l- SA{t)) + fdx'fAix', t) . (27) 

The above formula yields 

(28) 
In order to obtain a tagged particle survival probability 

S^.j^{t) we substitute expressions (l28l) into Eq. (|4]). In 

the long-time limit, the leading term of the resulting sum 

reads 



^^:^W-^(f)^A(t) 



(29) 



The asymptotic survival probability ([29)1 decays ex- 
ponentially with time, the decay rate being indepen- 
dent of n. The hard-core interaction manifests itself 
only through the n-dependent prefactor. Moreover, the 
asymptotic decay rate of Sf^.j^{t) is the same as that for 
the single noninteracting particle. This can be under- 
stood on physical grounds. Due to the mutual entropic 
repulsion among the particles, it is highly probable that 
the long-lived tagged particle eventually remains alone 
in the interval. Therefore its long-time dynamics should 
resemble that of the noninteracting particle. 



C. Fixed initial concentration of particles 

Consider now fixed c and random N case. The survival 
probability of the n-th particle, S:^{t), is obtained from 
S^.j^{t) by performing the summation (cf. Eq. ((1 



s:{t) 



N=n 



Stj,{t)P{N) 



(30) 



where P{N) is the Poisson distribution ([5]). In the long- 
time limit, all terms of the above sum are proportional to 
S'A(t) (cf. Eq. dH])). The asymptotic survival probability 
is given by 



s^{t) 



r (cX/2)"-l _ 

^^ (n-D! ' 



cL/2 



SAit) 



(31) 



The above result is formally similar with Eq. (f24|) . How- 
ever, the physical background of both asymptotics is 
rather different. In contrast to Eq. (|24p . Eq. ([5T|) rep- 
resents the overall contribution of all possible initial par- 
ticle numbers for which n < N (while Eq. ([M)) gives only 
the first term of the sum 



D. Random interval length 

Finally, let us derive the asymptotic tagged particle 
survival probability when L is drawn from the exponen- 
tial distribution with the mean 1/A. This corresponds 



to the model of diffusion on the real line with randomly 
distributed perfectly absorbing traps. The concentration 
of traps is uniform and equals A. We assume that the ini- 
tial concentration of particles c is also uniform. Actually, 
for non-interacting particles, there exists a large body of 
literature on the subject, see e.g. Refs. l24l. \2a. I33l - l37l and 
references therein. 

The average of the survival probability S:^{t) over the 
probability distribution of the interval lengths is given by 



s^it) 



dLXe-'^'^S^it) . 



(32) 



We substitute the asymptotics (|3ip into the above inte- 
gral. After that the integral is evaluated by a saddle- 
point method (see Appendix IBJ) . The resulting long-time 
asymptotics of S^ (i) reads 



S^it) ~ CiDt)'^^''+^^/^ejip< 



where the prefactor is 



27r2(A-f-) Dt 



1/3^ 



C = 



27rl6A (c/2)" (^2TT^){2n+l)/6 



(n-l)!(A-Kc/2)("+2)/3 ■ 



(33) 



(34) 



As it is for the single noninteracting particle;^ the 
decay of the survival probability is slower than its de- 
cay when the interval length is fixed. The stretched ex- 
ponential relaxation in Eq. (j33l) results from very large 
(and extremely rare) interval length L which enhance the 
asymptotic long-time survival probability.'^* A remark- 
able feature, which is not presented in a noninteracting 
model, is the dependence of the factor in the exponential 
on c. That is on the concentration of particles and not 
only on the concentration of traps A. This dependence 
is the consequence of the entropic repulsion among par- 
ticles - the higher the initial concentration the stronger 
the effective force that pushes the tagged particle into 
the trap. 



V. CONCLUDING REMARKS 

The main results of the paper are the asymptotics (|2ip , 
(|M)1 . dH]), dnH), dMl)- When the right boundary is reflect- 
ing, the particles located between the reflecting bound- 
ary and the n-th tagged particle create an effective (en- 
tropic) repulsive force which shortens the time spent by 
the tagged particle in the interval. This is reflected in 
the dependence of the decay rate of asymptotic tagged- 
particle survival probability (|2T]) on iV — n. If the initial 
concentration is given instead of the precise initial num- 
ber of particles, then Eq. (IMl) yields the asymptotic sur- 
vival probability averaged over all possible initial particle 
numbers. In fact, Eq. (p4|) represents only the first term 
of the sum (P5)l . Thus the long-time properties are con- 
trolled by the initial condition in which the n-th particle 
is the right-most one. 
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The situation is rather different when both boundaries 
are absorbing. For a given N, the long-time regime is 
governed by the situations when the tagged particle re- 
mains alone in the interval. Then its dynamics resembles 
that of the noninteracting particle. The only evidence of 
the hard-core interaction is the pre-exponential factor of 
the asymptotic survival probability (j29p . For a given c, 
all terms of the sum (j30p are asymptotically proportional 
to the survival probability of the single noninteracting 
particle S^it). Hence all initial conditions with all possi- 
ble initial particle numbers contribute to the asymptotics 
(Uni)- When L is random (and c is given), the tagged- 
particle asymptotic survival probability p3p exhibits a 
stretched-exponential decay. This is in parallel to what is 
known for the case without interaction.^^ A remarkable 
consequence of the hard-core interaction is the depen- 
dence of the asymptotic survival probability ((33)) on the 
combination X + c/2. That is, not only the concentration 
of traps A shortens the lifetime of the tagged particle but 
also the initial particle concentration c. 

Finally, notice that the results obtained in the present 
paper (and in Ref. \3W are different as compared to the 
result obtained in Ref. [23. The difference stems from 
the fact that in the present study (and in Ref. [Sy) all 
particles are absorbed at the boundary whereas in Ref. 
[271 only the tracer particle "senses" the presence of the 
absorbing boundary. 



dF{x,t)/dx, the integration over x is carried out which 
resuhs in Eq. ([T^. 



Appendix B: Laplace's approximation 

To keep the paper self-contained let us present an 
asymptotics of the integral arising in Sec. lIVDf S. 



I{t) = I dL L" exp 
/o 



-aL 



f3t 
I2 



a,/3>0. (Bl) 



The substitution L — t^^^x brings us to the expression 



f-OO 

lit) = <("+i)/3 / dx x" exp 
Jo 



_^l/3 



P 



t''nax + — 
x^ 



(B2) 

The above integral is approximated by the Laplace's 
method.— The approximation of the integral reads 



lit) ^ C„ i(2"+l)/6 



exp 



-- {2l3ah 



24\l/3 



where 



L'n — 



27r (2/3)(^"+i)/6 
3 a("+2)/3 



(B3) 



(B4) 
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Appendix A: Derivation of S'„:jv(t) and S„{t) 

Notice that f{x,t) — dF{x,t)/dx. By definition we 
have (for the sake of brevity the arguments of F{x, t) are 
omitted) 



Sn:N{t) — 



N' 



N- 



dF 



dxi^"-Ml-F)""":— . (Al) 
ox 



Next we expand the term [1 — (1 — F)]" according to 
the binomial theorem. We get 



ri-l 



^«:.(t)=nnE(-i)'^^"'' 



^ ' fe=0 
(-1) 



k 



(A2) 



f^ d 



N-n+k+l 



which, after performing the required integration, yields 
Eq. Q. 

Derivation of the survival probability (|12|) requires sim- 
ilar steps. First, we expand the exponential in Eq. ([7]) 
into the power series, then, using the fact that f{x,t) — 
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